Antiphase Synchronization in a Flagellar-Dominance Mutant of Chlamydomonas 
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Groups of beating flagella or cilia often synchronize so that neighboring filaments have identical 
frequencies and phases. A prime example is provided by the unicellular biflagellate Chlamydomonas 
reinhardtii, which typically displays synchronous in-phase beating in a low- Reynolds number version 
of breast stroke swimming. We report here the discovery that ptxl, a flagellar dominance mutant of 
C reinhardtii, can exhibit synchronization in precise antiphase, as in the freestyle swimming stroke. 
Long-duration high-speed imaging shows that ptxl flagella switch stochastically between in-phase 
and antiphase states, and that the latter has a distinct waveform and significantly higher frequency, 
both of which are strikingly similar to those found during phase slips that stochastically interrupt 
in-phase beating of the wild type. Possible mechanisms underlying these observations are discussed. 

PACS numbers: 87.16. Qp, 87.18.Tt, 47.63.-b, 05.45.Xt 



Living creatures capable of motion seldom restrict 
themselves to a single mode of self-propulsion. Pairs 
of appendages of multilegged organisms can be actuated 
synchronously in-phase, synchronously out-of-phase, or 
asynchronously, typically under neuronal control through 
a "central pattern generator" [1 . In the world of aquatic 
microorganisms, where there is no direct analog of a cen- 
tral nervous system, the cilia and flagella adorning algae 
and bacteria are the "limbs" which exhibit various sy- 
chronization modes, generating swimming [2 . Within 
a given eukaryotic organism, the undulations of flagella 
which arise from molecular motors distributed along the 
length of the filaments can be found to synchronize in 
two stereotypical ways. Biflagellate cells epitomized by 
the alga Chlamydomonas [3 display synchronous beating 
with identical frequencies and phases [H [5] . Those with 
multitudes of cilia or flagella, such as unicellular Parame- 
cium [6 or multicellular Volvox [7 , exhibit metachronal 
waves in which flagella with a common frequency have 
phases that vary monotonically with position. Theory 
[8-10 suggests that these modes of synchronization can 
arise from fluid dynamical coupling between flagella, pos- 
sibly assisted by waveform compliance. 

Flagellar synchronization is more complex than the 
simplest deterministic models of coupled oscillators 
would suggest; beating is intrinsically stochastic, cells 
can switch between synchrony and asynchrony [5 , and 
flagella existing within a single organism can be func- 
tionally distinct. These features are well-established for 
Chlamydomonas; the flagella of wild- type {wt) cells typi- 
cally exhibit a noisy in-phase (IP) breaststroke (Fig. flk). 
Termed cis and trans for their proximity to the cell's eye- 
spot, the two flagella are differentially affected by internal 
calcium levels, exhibiting a tunable flagellar dominance 
pT] that allows for phototactic turning. 

We report here an alternative mode of synchronization 
not previously quantified |12j in eukaryotes, in which flag- 
ella lock in antiphase (AP) synchronization. For a range 



of experimental conditions [TT, this behavior can be sus- 
tained in time by the "flagellar-dominance" mutant ptxl 
of C. reinhardtii [H]. While vegetative cells of ptxl ex- 
hibit no gross motility defects, they have defective pho- 
totaxis [I2j [HI [15] thought to arise from lack of Ca^+- 
dependent flagellar dominance. We discuss mechanisms 
proposed for AP synchronization [HI [TBHTO] , and suggest 
that our observations support active filament models |20j 
which exhibit discrete undulating modes of flagella. 

Wild-type (GG125) and ptxl mutant (CG2894) strains 
of Chlamydomonas reinhardtii [21 were grown photo- 
autotrophically in Tris-Minimal medium [22 with revised 
trace elements [23 and air bubbling in a diurnal growth 




ptxl IV 



FIG. 1. (color online). Waveforms of C reinhardtii. 
Logarithmically-scaled residence time plots averaged over 
(9(10^) beats overlaid by waveforms during a cycle, color- 
coded in time. The wild- type [wt) displays in-phase (IP) 
breaststroke beating (a) stochastically interrupted by phase 
slips (b) in which one flagellum (here, trans) beats faster with 
an attenuated waveform. The mutant ptxl displays an IP 
state (c) nearly identical to the wild- type (a), and a high- 
frequency antiphase (AP) state (d). Large and small ovals 
indicate cell body and eyespot, respectively. 
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FIG. 2. (color online). Beating dynamics, (a) Phase difference A = {iptrans — '0cis)/27r showing half-integer jumps between IP 
and AP states. Insets show waveforms in the two states, (b) Instantaneous frequencies of AP and IP states, (c) Distribution 
of per-beat instantaneous frequencies during IP beating and of the faster flagellum during slips, across all sampled wt cells, (d) 
IP and AP per-beat instantaneous frequency distributions, across all sampled ptxl cells. 



chamber at 24° C on a 14:10 h light-dark cycle with a light 
intensity of 90 /iEm~^ s~^ [5j. Cells were harvested from 
1 or 2 day-old cultures at a density ~ 6 x 10^ cells/ml, 
during hours 4 and 5 of the subjective day. Cells were 
washed in fresh buffer HKC-40/10 (5 mM HEPES, 40 
mM KCl, 10 mM CaCl2, pH 7.2) and allowed to regrow 
flagella for at least 2 hours. Density and motility were 
monitored prior to harvesting and after washing. Spe- 
cially designed cylindrical PDMS chambers 15 mm in di- 
ameter x4 mm height were cast in custom aluminum 
molds and plasma-etched onto 22 x 50 mm cover slips. 
Chambers were placed on a Nikon TE2000-U inverted 
microscope with a x63 Plan-Apo water-immersion ob- 
jective (Carl Zeiss AG, Germany). Micropipettes were 
used to hold and orient cells as described previously [5|. 
Bright-field illumination was carried out using a halo- 
gen lamp with a red long-pass filter (> 620 nm) to min- 
imize phototactic behavior during experiments, which 
were performed in the absence of background illumina- 
tion. Video microscopy was performed at 1, 000 fps (Fast- 
cam SA3, Photron, USA), post-processed with custom 
MATLAB code. After each recording the filter was re- 
moved to locate the orange-colored eyespot and thereby 
identify the cis and trans fiagella. Based on the experi- 
ments with wt cells we concluded Chlamydomonas need 
to be acclimated for at least 20 — 30 min before char- 
acteristic synchronized breaststrokes are observed [H |5] . 
Data from 10 wt cells and 12 ptxl cells were analyzed. 

There are four key experimental results. The first is 
the existence of the AP state itself (Fig. [l]i), visualized 
here by discrete waveforms within one cycle, color-coded 
in time and overlaid on a spatial map of flagellar resi- 
dence time, averaged over many cycles. For reference, 
compare this to Fig.[l^, which shows the wt IP breastroke 
waveform. Here, the flagella simultaneously execute ex- 



tended "power strokes" followed by high-curvature "re- 
covery strokes", in which they are drawn forwards with 
distal portions sliding past the cell body. In the AP of 
the mutant, distinct power and recovery strokes are still 
clearly discernible, but as one flagellum executes the for- 
mer, the other proceeds through the latter. The mu- 
tant also displays an IP state (Fig. Ilfc) that is nearly [13] 
identical to the wt IP. For example, the areas A^p'-^ ^ 
swept out by the flagellum in both cases (i.e. the areas 
within residence-time plots in Fig. [l]) agree to within 
1%. In the case of ptxl^ evident also is the drastic re- 
duction in spatial extent spanned by both flagella during 
AP relative to the wt IP mode. This alteration of beating 
waveform occurs concomitantly with an abrupt increase 
in beating frequency, which together comprise our sec- 
ond observation. We extract flagellar phases ^l^ds, trans 
from Poincare sectioning of the dynamics as done previ- 
ously [5] , and define the interfiagellar phase difference as 
A = {iptrans " '0cis)/27r. For a typical ptxl cell. Fig. [2^ 
tracks A(t) over ~ 40 s. We see that A fluctuates around 
half-integer values during AP, but around integer values 
during IP. As seen in Fig. |2J our third finding is pre- 
cisely that flagella of ptxl stochastically transition back 
and forth between these IP and AP modes, in a manner 
reminiscent of the synchronous/asynchronous transitions 
of the wt [5] . Fig. [2J3 shows that the instantaneous beat 
frequency is indeed higher in AP (z^ap • 82 ±4 Hz) than in 
IP {uip : 58 ± 5 Hz). Fourth and finally, we highlight the 
striking similarities between the AP state and the state 
of the flagellum that accumulates one additional cycle 
during a phase slip of the wt [5]. This is evidenced by 
the equivalence both of the waveforms (Fig. [l]3,d, areas 
^siL'^^AP ^gi"66 to within 5%), and of the frequencies 
(Fig. [2fc,d). The latter figure showing also that wt and 
ptxl cells beat at similar frequencies during IP. 
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FIG. 3. (color online). Synchronization dynamics. Phase 
plane of polar angles Ods, trans of a single point on each flag- 
ellum reveals the IP (green) and AP (red) synchronization of 
a single cell (a), and (b) the average over 6 cells, averaged 
over (9(10^) beats and resampled at 15 points, equally-spaced 
in time. Shaded regions in (a) indicate one standard devia- 
tion of stochastic fluctuations. (c),(d) are sample timeseries 
for evolution of Ods, trans during a typical transition event, (e) 
Phase difference dynamics during transitions from AP to IP 
(orange) and IP to AP (blue) for 60 transitions taken from 
all sampled cells, with means (solid lines) and standard devi- 
ations (shaded). Transitions have been vertically aligned by 
plotting difference modulo 1. Dashed lines are fits to data. 



The hypothesis that there is a second, distinct beat- 
ing mode of the flagellum is supported by estimates of 
the flagellar force F and power P ^. In a caricature 
of the power stroke we imagine a straight flagellum of 
length L pivoting from initial polar angle ^o to a final 
one Of during half the beat period. Using resistive force 
theory we integrate the normal component of the viscous 
force along the filament to obtain F ~ 2(±uA^ where 
C_L is the perpendicular drag coefficient per unit length 
and A is the waveform area defined previously. A sim- 
ilar calculation yields the power P ^ (2/3)FV, where 
V = LO is the tip speed of the flagellum. Ratios of the 



product uA thus serve as measures of relative force in 
different beats. Restricting to a subset of cells whose 
flagella were most planar, averaged values of the pairs 
(z/. A) for the four states of interest are: ptxl IP: (57.2 Hz, 
147.3 /im^), ptxl AP: (81.0 Hz, 105.1 /im^), wtIP: (59.4 
Hz, 148.8 /im^), wt slip: (82.0 Hz, 110.1 /im^). We find 

^ip'^V^Ip^ = 0-99 ± 0.06 and Fj^V^.Tip = 0-98 ± 0.07. 
The quantitative match of these ratios supports the iden- 
tification of a wt slip with the transient appearance of a 
higher mode, and the fact that the common value is ac- 
curately unity would also imply equal force generation in 
the two states. Intriguingly, the ratio of the average AP 
and IP frequencies for ptxl and of the average slip and 
IP frequencies of the wt are nearly identical, with a value 
close to 4/3. Finally, detailed studies [24] show that the 
peak force during IP power strokes are ^ 20 pN with 
peak powers (9(10 fW). These are in agreement with es- 
timates from time-resolved PIV measurements of energy 
dissipation in the fluid around free-swimming cells [25]. 
The polar angles {0 ds •, trans) nieasured from the cell 
midline to equivalent points on the two flagella define 
a low-dimensional phase space with which to quantify 
synchrony. Figures [3^,b show IP and AP motion in this 
space for a single cell and a multi-cell average. Individual 
cells orbit fairly close to the diagonals, but the mean dis- 
plays remarkably precise IP and AP motion, with phase 
coherence maintained during power and recovery strokes. 
Transitions to and from these two types of synchrony 
(Figs. |3]3,d) are always initiated by one flagellum, either 
cis or trans^ which undergoes alteration of beating mode 
first [T3 . Using Poincare sections we examine the re- 
emergence of synchrony during transitions between the 
two modes using the difference (V^iead — '0foiiow)/27r be- 
tween the phase of the flagellum that leads the transition 
and that which follows. The transition dynamics of re- 
spectively AP^IP and IP^AP obey an equivalent func- 
tional form derived, on a phenomenological level, from a 
noisy Adler equation for A [5 



K = -v\K)+m- 



(1) 



Here V{/S) = —SuA + ^(A), with 6u an intrinsic fre- 
quency difference and U an effective potential periodic in 
A, and ^{t) is a noise term. Applying this to either type 
of synchrony in ptxl we expect 5iy c:^ due to the lack of 
flagellar dominance [T^. The most parsimonious model 
would then he U = -ecos(27rA), with e > for AP^IP 
and e < for IP^AP. Solving for the deterministic dy- 
namics (^ = 0) in a scaled time s = v{t — ti) centered at 
the inflection point of the transition t^, where v is the IP 
frequency, we obtain A = — (l/27r) cos~^ tanh(5/r), with 
rescaled relaxation time r = l/(47r^e/i/). Fits to the data 
yield tap^ip = 1.65 ± 0.02 and tip^ap = -2.07 ± 0.04 
(Fig. [3]3,d) and thus cap^ip/^ — 0.015 and cip^ap/^ — 
—0.012, consistent with the wt [5 . 

The necessity to invoke couplings of opposite sign to 
account for the AP and IP states within the simplest 
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FIG. 4. (color online). Synchronization mechanisms. In the 
elastohydrodynamic model, co-rotating spheres (a) synchro- 
nize in phase, while counter-rotating spheres (b) do so in an- 
tiphase. Top row: motion of sphere 1 at two possible phases 
01 induces flows (blue arrows) which alter the trajectory of 
sphere 2, either speeding it up (green), or slowing it down 
(red). Bottom row: the converse perspective, (c) Elastic cou- 
pling between flagella can induce either IP or AP modes. 



Adler equation ([T]) provides a natural starting point for 
a discussion of mechanisms proposed for synchroniza- 
tion. Two key issues arise: the structure of the poten- 
tial U and the origin of the coupling constants. With 
^z/ = 0, the solution to the Fokker-Planck equation 
for the probability distribution function P(A) associ- 
ated with gives /3U = -log[P(A)] with /3 related 
to the noise in the usual manner. The function PU 
so determined [26] will be a bistable potential with lo- 
cal minima at integers and half-integers. This could be 
accommodated by higher-order Fourier components, as 
U{A) c^ -ecos(27rA) - acos(47rA), with e > and 
a > e/4. An alternative to this picture of a fixed potential 
landscape U{A) with stochastic hopping between locally- 
stable minima is a fluctuating landscape switching be- 
tween potentials Ujp and Uap^ the former with minima 
only at integers, the latter at half-integers. Within the 
limitations of a phase-oscillator description in which am- 
plitude dynamics are suppressed, the distinction between 
these views is fundamentally a matter of which degrees 
of freedom are considered part of the dynamical system 
and the relative time scales for those variables. In fact, 
precedent for a fluctuating landscape can even be seen in 
the 1^;^ [5] , in which asynchronous beating ( "drifts" ) cor- 
responds to a washboard potential tilted by a large Su 
so there are no local minima, while synchronous beating 
has 6u small enough to allow local minima. 

Models of synchronization based on hydrodynamic 
coupling often represent flagella by microspheres exe- 
cuting trajectories driven by a tangential internal force. 
That force may be considered constant along a trajectory 
with elastic compliance [9|, or the trajectories are rigid 
and the forcing varies with orbital phase [8]. The mech- 
anism of synchronization in the first class is illustrated 
in Fig. [4|i,b. Measuring the phase angles (0i, ^2) of the 
spheres as indicated, cilia would be modelled by orbits 
in the same direction, say cji = ^i > and uj2 = (J)2 > ^ 
(Fig. pa). If sphere 1 lags 2 then the flow produced by 



1 will push 2 to a larger radius. If the internal force is 
constant, ^2 will decrease, and 1 will catch up. Con- 
versely, if 1 leads 2 then it pushes 2 inward, so 2 acquires 
a higher phase velocity and will catch up. Similarly, the 
flow induced at 1 by 2 leads to consistent results, showing 
that co-rotating IP motion is stable. To model Chlamy- 
domonas the two spheres must be counter-rotating, with 
say uji > and UJ2 < (Fig. [4J3). Then, these considera- 
tions, together with anisotropy of the stokeslets, predict 
stable AP synchronization. Indeed, the coupling con- 
stant in (IT]) scales as e oc —ujiuj2 and is negative (positive) 
for CO- (counter-) rotation. In this simple model the AP 
beating of ptxl is the 'normal' behavior, and the IP mode 
is anomalous! The situation is not so clear, though, for 
if the relationship between orbital radius and phase ve- 
locity is reversed then the coupling changes sign J6l [17] . 
This relationship could be influenced by mechanosensi- 
tive cues [27]. In the class of models with forcing that 
varies with phase angle, synchronization can be under- 
stood by similar means in terms of the flow induced by 
one sphere at the other. Allowing for non-circular tra- 
jectories as well as proximity to a no-slip surface leads to 
the possibility of an effective potential with the higher- 
harmonic structure discussed above, stabilizing both IP 
and AP patterns [81 [19] . The difficulty in determining 
the relevance of these arguments to ptxl is precisely that 
the two modes of synchronization are associated with 
very distinct waveforms, with potentially different com- 
pliances, internal forcing, and proximity to the cell sur- 
face. A third model [18 builds on the fact that transient 
deviations from locked phases will lead to yawing mo- 
tion of the cell body which can produce differential forces 
on the flagella, bringing them back into phase. While 
such a mechanism may pertain to free-swimming cells, it 
is not immediately clear how it can encompass the ap- 
pearance of both IP and AP states of cells held strongly 
on micropipettes, where we observe only minute angu- 
lar displacements (below 1° in both states). The pres- 
ence of the cell-body itself appears not to be essential for 
synchrony of the two flagella, for a wt-\ike breaststroke 
gait has been observed in isolated flagellar apparati (ax- 
onemes still connected through their basal bodies), after 
reactivation by ATP [28] . 

No existing models of eukaryotic flagella explain the 
antiphase waveform nor explain why it should emerge. 
Approaches based on optimizing swimming efficiency or 
nutrient uptake in a model of Chlamydomonas [29 do 
find a mode comparable to the IP state. Perhaps the 
AP waveform is not optimal in any conventional sense, 
but instead exists as one of a discrete number of modes 
that can emerge from sliding filament models [20 . It will 
be important to establish whether the higher frequency 
and distinct waveform are properties intrinsic to a single 
flagellum, or derive from interactions between the two; 
key insight will likely be gained from examining flagellar 
dynamics of uniflagellated double mutants of ptxl. 



The physiology of stochastic transitions in the pattern 
of flagehar beating (i.e. shps or transitions in and out of 
AP) is currently unknown; we hypothesize that fluctu- 
ations in the concentration of a small molecule or ion 
might be the origin. One candidate would be Ca^+, 
which in isolated and reactivated flagellar axonemes is 
known to control the waveform [30]. Interestingly, cal- 
cium ions are also responsible for the contractility of stri- 
ated fibers that connect the basal bodies of the two flag- 
ella [31] , which in turn may act as a spring with variable 
stiffness. The current state of this potential spring may 
influence the preferred mode of synchronization. Indeed, 
generalizing the orbiting-sphere model [9^ to include an 
elastic connection between flagella bases can lead to sta- 
bilization of either the IP or AP modes (Fig. pfc), de- 
pending on microscopic details of the elasticity. In the 
simplest linear spring, for example, the AP mode (termed 
'parallel' by Riiffer and Nultsch [12j) can be selected, 
for it is the mode in which the relative displacements of 
the flagellar connections within the cell body are most 
nearly constant. The role of these fibers for flagellar syn- 
chronization may be clarified by altering their mechanical 
properties by chemical or other means. 
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